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Abstract 

Oh! 

r-| ! A minimal supersymmetric SO(IO) model with one 10 and one 126 Higgs superfield has 

recently been shown to predict all neutrino mixings as well as the solar mass difference squared 
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in agreement with observations. Two assumptions critical to the predictivity and success of 
the model are that: (i) the superpotential includes only renormalizable terms, thereby limiting 
the number of free parameters and (ii) the triplet term in the type II seesaw formula for 
neutrino mass dominates, leading to the sum rule = c{Mii — Mg) that is responsible for 
large mixings. However, CKM CP phase is constrained to be in the second or third quadrant 
requiring a significant non-CKM component to CP violation to explain observations. We revisit 
this issue using type I seesaw formula for neutrino masses and obtain the following results: (i) 
we show that the above sumrule responsible for large mixing angles can also emerge in type I 
seesaw models; the detailed predictions are however not compatible with present data for any 
choice of CP phases, (ii) We then show that addition of a nonrenormalizable term restores 
compatibility with neutrino data and CKM CP violation both in type I and type II cases. We 
further find that (iii) the MSSM parameter tan/3 > 30 in the type I model and (iv) lepton flavor 
violation and lepton electric dipole moments which are accessible to proposed experiments in 
both type I and type II models. We also discuss the unification of the gauge couplings in type 
I model which requires an intermediate scale. 



1 Introduction 



The observations involving the solar and atmospheric neutrinos together with those using the 
accelerator and reactor neutrinos have now conclusively established that neutrinos have mass 
and they mix among themselves. In conjunction with the negative results from CHOOZ and 
PALO- VERDE reactor experiments, a reasonably clear outline of the mixing pattern among 
the three generations of neutrinos has emerged. Of the three angles needed to characterize 
these mixings, ^12, ^23 and ^13, the first two, responsible for solar and atmospheric neutrino 
deficits respectively, are large, and the third corresponding to reactor neutrinos is small 
One of the major experimental issues in neutrino physics now is to make the knowledge of these 
angles more precise. 

Foremost among the theoretical challenges that have already emerged from these discoveries 
are, first an understanding of the smallness of neutrino masses and second, understanding the 
vastly different pattern of mixings among neutrinos from the quarks. Specifically, a key question 
is whether it is possible to reconcile the large neutrino mixings with small quark mixings in 
grand unified frameworks that unify quarks and leptons. 

The first challenge, i.e. the lightness of neutrino masses is elegantly answered by the seesaw 
mechanism |2] which requires an extension of the standard model that includes heavy right 
handed neutrinos. The light neutrino mass matrix is obtained by integrating out heavy right 
handed neutrinos and one gets 



where Mj; is the Dirac neutrino mass matrix and M/j is the right handed Majorana mass 
matrix. It is therefore important to explore whether one can answer the second puzzle of large 
neutrino mixings within the seesaw framework. 

The above formula for the neutrino mass matrix is called the type I seesaw formula. The 
right handed Majorana mass scale, Vn-, is almost determined by the mass squared difference 
needed to understand the atmospheric neutrino data to be around 10^^ GeV, if we assume 
that the Dirac neutrino mass is same as up-type quark mass. Before proceeding to discuss the 
implications of this large right handed neutrino mass, let us discuss the nature of the seesaw 
formula. 

We will consider a class of models where the right handed neutrino mass is not put in by 
hand but arises from a renormalizable coupling of the form JMNAr, where is a right handed 
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neutrino, / is a coupling constant and A/j is a Higgs field whose vacuum expectation value (vev) 
gives mass to the right handed neutrino. This is a natural feature of models with asymptotic 
parity conservation, such as those based on SU{2)i x SU{2)r x U{1)b~l or any higher gauge 
group such as SO(IO)), where the field is part of an SU{2)ii triplet field. Parity invariance 
then implies that we also have an Ji/uAl coupling term as a parity partner of the NNAn 
coupling. In this class of theories, whenever A/j acquires a vev, so does A^ and they are related 

2 

by the formula (A/,) = = where is a weak scale and is the A^ vev and A is a 

coupling constant in the Higgs potential. The A^ vev contributes a separate seesaw suppressed 
Majorana mass to the neutrino leading to a modified seesaw formula for neutrino masses given 
below. 

M'^ = M,-M,^M^\M^r, (2) 

where = /f l and M/j = fvR. This formula for the neutrino mass matrix is called the type 
II seesaw formula j3]. In the case where right handed Majorana masses are heavy enough, the 
second term in the type II seesaw formula can be negligible, and the first term, = fv^, 
is dominant. We will call this pure type II seesaw. When both terms are comparable, we will 
call this mixed type II seesaw. 

Coming to the large scale for right handed neutrino mass (i.e. 10^^ GeV or higher), we note 
that it suggests supersymmetric grand unified theory (GUT) as a natural framework to study 
neutrino masses, since experimental data suggests the gauge coupling unification scale in the 
minimal supersymmetric extension of the Standard Model (MSSM) to be, Mq ~ 2 x 10^^ GeV, 
which is close to the seesaw scale. The minimal grand unification model for neutrinos is the one 
which is based on the SO (10) group since all standard model fermions and the right handed 
neutrino fit into the 16-dimensional representation of SO(IO), resulting not only in a complete 
unification of the quarks and leptons but also yielding possible relations between the quark and 
lepton mass matrices. One may therefore hope that the neutrino oscillation parameters might 
be predictable in an SO(IO) theory. 

There are two simple routes to realistic SO(IO) model building. In the first class, one may 
have smaller representations for the Higgs fields like 10 and 16 multiplets. In this case, the 
nonrenormalizable terms are added to the superpotential to implement the seesaw mechanism. 
These models have the disadvantage that they break R-parity which then induces rapid proton 
decay at an unacceptable level. 

An alternative is to introduce both 10 and 126 Higgs multiplets to give fermion masses. 
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In this class of models, the R-parity is an automatic symmetry of the model. This naturally 
prevents the baryon and lepton number violating terms that give rise to rapid proton decay 
and also guarantees a naturally stable supersymmetric dark matter. 

In SO(IO) models of this type, both the and Ar fields are part of the 126 multiplet. 
The above mentioned couplings that contribute to seesaw formula arise from the couplings of 
16 matter spinors to 126 multiplet. If the rest of the Higgs sector is appropriately chosen, this 
leads to type II seesaw formula for neutrinos. 

An interesting class of renormalizable SO(IO) models with 126 was proposed in Ref.jS] 
where it was shown that if in addition to the Higgs multiplets that break SO(IO) and do not 
couple to fermions, one chooses only one 10 and one 126 Higgs field, then the model provides a 
completely realistic description of the charged fermion sector of the standard model and is very 
predictive in the neutrino sector without any extra symmetry assumptions. The reason for this 
is that all the Yukawa parameters of the model are determined by the quark and charged lepton 
sector. This model is called the minimal renormalizable SO(IO) model, since unlike the SO(IO) 
models with 16 Higgs fields, it does not add any nonrenormalizable terms to the superpotential. 
The presence of the 126 representation that allows the neutrino flavor structure to be related 
to other fermion mass matrices is at the heart of the predictivity of the model. 

In the Refs. jnj[7||H], this model was analyzed using the type I seesaw formula to see whether 
the neutrino oscillation parameters predicted would agree with observations. They found that 
the atmospheric and solar mixings can be large, and the phenomenological predictions are 
studied jH] . However they found that mass squared ratio for the solar and atmospheric data is 
Amg/Am^ ^ 0.2. This result is incompatible with the recent combined data analysis of the 
solar mixing angle and mass squared ratio [2]. 

A new approach to discussing neutrinos in this model was presented during the past year. 
Considering only the second and third generation sector, the authors of Ref. [TU] pointed out 
that the type II seesaw in the minimal SO(IO) model can provide a natural way to understand 
the maximal atmospheric mixing due to the convergence of bottom quark and tau lepton masses 
when extrapolated to GUT scale. The reason for this that in the the pure-type II seesaw case, 
the 126 Higgs coupling leads the neutrino mass matrix to be proportional to — Mg, where 
Md is a down-type quark mass matrix and Me is a charged-lepton mass matrix. It is then easy 
to see that h-r mass convergence makes the 3-3 entry of the neutrino mass matrix very small, 
leading to large atmospheric mixing angle. While this observation was interesting, a priori, it 
was not clear if this would survive once the model is extended to include three generations. It 
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was however shown in Ref . ^T] , that the same b-r mass convergence that led to large atmospheric 
mixing, also leads to large solar mixing while keeping the 1-3 mixing angle small as required by 
data. It also predicts the ratio of solar to atmospheric mass squared difference i.e. Amg/Am^ 
to be in the right range. This establishes that the minimal renormalizable SO(IO) models with 
126 Higgs provide a very interesting way to understand neutrino masses and mixings within a 
complete quark-lepton unified framework. 

In Ref.|llj. all the CP phases, including the Kobayashi-Maskawa (KM) phase, were set 
to zero (or 180 degree) and it was assumed that all CP violating effects owe their origin to 
the SUSY breaking sector. In Ref. [12], the pure type II model including all CP phases was 
analyzed, and it was found that one can maintain the predictivity of the model despite the 
appearance of the phases; but the KM phase must be in the second quadrant in the p-r] plane 
to maintain the neutrino predictions (where we have adopted the Wolfenstein parameterization 
for the CKM matrix). This implies that the model is substantially different in CP violating 
sector from the standard CKM model jT2] where the KM phase must be in the first quadrant. 
Such non-CKM CP violation, for instance, could be in the squark mixings. Phenomenologically 
speaking, there is nothing wrong with this possibility although admittedly it will require many 
random parameters to reproduce the observed data on CP violation. 

Since the CKM model is simple and seems consistent with known data on CP violation and 
it would be interesting to see if the above predictive model for neutrinos can coexist with the 
simple CKM CP violation. One of the objectives of this paper is to study this using the type 
I seesaw formula. It was mentioned in Ref.p2] that the main reason for the constraint on CP 
phase being in the third quadrant in the pure type II model has to do with fitting the electron 
mass and if a nonrenormalizable term is introduced to fix this problem, then the CP phase 
could be of CKM type. 

In this paper, we focus primarily on the same minimal SO(IO) model but with both type 

I and type II seesaw formula for neutrinos. We obtain the following results: (i) We first show 
that in a certain limit, one can get the relation ~ c{Md — Me), so that the basic mechanism 
that led to the maximality of solar and atmospheric neutrino mixings in the case of pure type 

II seesaw formula is preserved in the type I case. As far as we know, this fact was not noticed 
in any of the earlier papers that analyzed type I seesaw in minimal SO(IO). We then show 
that the model is in conflict with neutrino data with or without CP violating phases. This 
is different from the pure type II case where, one could get the neutrino predictions correctly 
within the renormalizable framework without any CP phase in the Yukawa couplings as well 
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as the KM phase in second quadrant, (ii) Secondly, we show that we can restore compatibihty 
with KM CP violation for both type I and type II cases by using the SO (10) model as an 
effective theory at the GUT scale, where we generate a specific nonrenormalizable term. This 
not only leads to successful predictions for neutrino mixings but also maintains the CKM model 
for CP violation. Of course, the resulting model is not a minimal model anymore since at the 
Planck scale, it can emerge from a renormalizable theory with two pairs of 126 fields (rather 
than one in the minimal case), (iii) We also find that the type I model works only for values 
of tan/3 higher than 30, providing a way to test the model in supersymmetry experiments, (iv) 
We evaluate the lepton flavor violating processes and electron electric dipole moment (EDM) 
by using the predicted leptonic mass matrices in both type I and pure type II case and show 
that they are in a range accessible to the present experiments, (v) We find that for the type 
I model to work, intermediate scale (vr) must be below the GUT scale. We present a brief 
discussion of unification of gauge couplings in the model to show that this can indeed happen. 

This paper is organized as follows: in section 2, we give a basic framework of the minimal 
SO(IO) model and show that the relation M^, ~ c{Ma — Mg) needed to understand large 
neutrino mixings follows in a certain limit even with the type I or mixed type II seesaw. We 
examine its predictions for the case where the type I seesaw formula is used for discussing 
neutrino masses and KM phase is kept around 60-70 degrees as required in the standard model 
fit. We find that the predictions for solar neutrino oscillation do not agree with the current 
combined data analysis of mixing angles and ratio of mass squared differences for this case. 
We also discuss the restriction on tan/3 in the mixed type II and type I models. In section 
3, to remedy the situation of neutrino masses and mixing angles, we include specific types of 
non-renormalizable terms in the superpotential and show how they are helpful in reconciling 
the neutrino predictions with CKM CP violation for both type I and type II models. We also 
discuss the origins of these terms. In sec. 4, we discuss the predictions of the model giving 
numerical results for various parameters as well as neutrino masses and mixings. In section 5, 
we discuss the gauge coupling unification in the model and show that the B — L scale required 
for obtaining neutrino masses is compatible with gauge coupling unification. In sec. 6, we 
evaluate the leptonic flavor violation and lepton electric dipole moment for both type I and 
type II models. Section 7 contains our conclusions. 
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2 Minimal Renormalizable SO (10) Model and Its Pre- 
dictions 



In the minimal SO(IO) model, Yukawa interactions are given as the couplings of the 16- 
dimensional matter spinor with only one 10- dimensional Higgs H and one 126 Higgs A. 
The superpotential for Yukawa interactions is written as 

Wy = ^hijtpitpjH + ^fijipiipjA . (3) 

The Yukawa couplings, h and /, are complex symmetric 3x3 matrices. 

The SO(IO) gauge symmetry breaks down by the Higgs mechanism. There are several 
breaking pattern in SO(IO) GUT. The SO(IO) symmetry is broken to the left-right group, 
^2231 = SU{2)l X SU{2)r X SU{3)c x U{1)b-l, using the 54 and 210 Higgses. To break 
SU{2)ji X U{1)b~l down to U{1)y, we will use 126 -|- 126 Higgs multiplets. 

The SO(10)-invariant superpotential, Eq.Q, includes MSSM Yukawa couplings plus right 
handed Majorana mass term which are written by using the MSSM superfields: 

Wy D Y.'^QiU^H^ + Yf^Q.D'^^Hd + Y^L.E^H, + Y'<L,N;rHu + ^^iVf iVJA^ , (4) 

where Hu and are MSSM Higgs doublets which are linear combinations of the SM doublets 
in 10 and 126 Higgs multiplets, and A/j is part of the 126 field A. We note that the VEV of Aji 
gives right handed neutrino masses and breaks SU{2)ji x U{1)b~l symmetry down to U{1)y- 
As noted the SO(IO) Higgs fields, H and A, contain two pairs of SU{2)l Higgs doublets, H^q 
and H^^, and the MSSM Higgs doublets are linear combinations of two pairs, 

= a.H^o + PdHt^,, (6) 

where la^.dP + I/^m,^^ = 1- The other linear combinations are assumed to be massive around 
GUT scale. The MSSM Yukawa couplings Y^''^''^''^ are given by linear combination of h and /, 
and the fermion mass matrices and Majorana mass matrix for right handed neutrino are given 

as"^ 

Mu = {h*au + ri3uW (7) 
= {h*aa + r(3d)va, (8) 
^The mass matrices are defined as —Cm — i^hMipR + ^[i^fi)^ Mnyjf + h.c. 
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where 



M, = {h*ad-3rPd)vd, (9) 
M,^ = - 3r/5„K, (10) 

Mn = rvR, (11) 

where Vu,d are VEVs of MSSM Higgs doublets and vr is a VEV of Ar. We denote that 
Vu = f sin/? and = f cos/3, where w = 174 GeV. Then we have sumrules for leptonic mass 
matrices: 

Me = CdMd + CuMu, (12) 

M,^ = ^^(M,-M,) + M., (13) 

Md- Me 

4cot/3 -4/3,//3d 

So if the quark mass matrices M„ and M^^ are our input, then leptonic mass matrices are 
determined by two complex parameters, c„ and barring the overall scale and the phase of 
Majorana neutrino mass matrix, Mr. Thus, using the masses of three charged leptons, we 
predict the neutrino mass matrices in terms of only one real parameter. This is an interesting 
feature of the minimal SO(IO) model and is responsible for the high predictive power of the 
model. 

Let us count the number of parameters in the fermion mass matrices in the minimal SO(IO) 
model. Rotating the fermion fields by a unitary matrix without violating SO(IO) symmetry 
(fermions in 16 are all rotating simultaneously), up-type quark mass matrix, M^, can be made 
real positive diagonal, M„ = diag ( m^j, mc, m^) = D^- Then down-type quark mass matrix, M^, 
is a general complex symmetric matrix, which has 6 complex parameters. The matrix can be 
written as M^ = UD^U^ where D^, is real positive diagonal, = diag(md, mg, rriiy), and [/ is a 
unitary matrix. The unitary matrix has 9 real parameters (3 angles and 6 phases) and we can 
parameterize the unitary matrix as f/ = P*UPd- The matrices P„ and Pd are diagonal phase 
matrices, P„ = diag(e*''^"/^, 6*'^'=/^, 1) and Pd = diag(e*'^''/^, e*''^^/^, e*'^''/^) and U is just same as 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix, Vckm, in which there are 3 mixing angles 
and 1 phase (KM phase). The phase matrices, Pu and Pd, are unphysical in the MSSM quark 
sector since they can be absorbed in right handed quark fields, but the phases are relevant 
parameters in the leptonic mass matrices through the sumrules in Eqs. (|12llT^ . One overall 
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phase of is irrelevant in the sumrules, using the rotation, — > e M^, M„ M^. This 
rotation is generated by {ad, Pd) (e^'^a^, e*'^/?^). So, we can choose the phase (pb to be zero. 
Finally, we have 19 parameters in the sumrules: 14 real parameters in quark mass matrices 
Mu and Md (6 quark masses, 3 mixing angles, 1 KM phase, and 4 unknown phases, (f)u,c,d,s), 2 
complex parameters, Cu and q, and 1 overall scale parameter in the Majorana mass matrix Mr 
(overall phase of M/j is not counted since it is not physical). On the other hand, the MSSM 
Yukawa couplings plus the right handed Majorana mass in Eq.(j3]), give us 32 parameters: 10 
parameters in and (6 quark masses and 3 mixings and 1 KM phase), 10 parameters in 
and in the same way, and then 12 parameters in the Majorana mass matrix which is 
generally a symmetric matrix. So, in the minimal SO(IO) model, we have less parameters than 
the MSSM. 

The quark masses, mixings and the KM phase are our inputs. We can redefine the basis 
of fermions, and the up- and down-type quark mass matrices are written as = D^Pl^, 
Md = VcKuDdPdycKM- The authors in Ref.[7] assumed that the phase matrices P^^ are real 
and they found a solution to generate the observed quark masses and the CKM matrix and the 
charged lepton masses in the case where (t)d,s,c = and = 0, vr. In our general analysis, we 
take those phases 4'd,s.u,c to be arbitrary. Those phases are constrained to obtain the bi-maximal 
neutrino oscillation data. 

Let us see how one can obtain a numerical fit of the parameters in the lepton sector at the 
GUT scale. The charged lepton mass matrix is written as = CdMd + c^M^. So, r and /i 
masses are approximately (neglecting generation mixing and phases), 

±mr ^ Cdmb + CulTlt, ±m^ CdTUs + CulTlc. (16) 

Using the relations rric/mt -C nig /mi, and m,- ~ m;,, we obtain Cd — -^m^/ms and c„ ~ 
~mh/mt{cd ± 1), and thus \cd\ is about 3, naively. In this naive approximation, however, the 
electron mass becomes about Cdfrid, which is too large, so we cannot neglect the off-diagonal 
elements, and we need a fine-tuning in det{Md + nMu), where k = Cu/cd, by a suitable choice 
of K. The numerical fit is given in Ref.jTj. The key ingredient to fit the electron mass is the 
strange quark mass. Actually, the current strange quark mass still has large experimental error 
and we can make it to be a parameter in the model. We show the values of the strange quark 
mass (at 1 GeV) and \cd\ by varying the KM phase in Fig.l. We take the KM phase to be 
in the first and second quadrant without loss of generality, since the sumrule does not change 
under the conjugation. The three real parameters are consumed to fix the three charged lepton 
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Figure 1: The strange quark mass (at 1 GeV) and \cd\ to fit charged lepton masses are shown. 

masses: n and \cd\ are determined and the phase of q is still undetermined. 

We now determine the neutrino mass matrix. First, in the pure-type II case, the light 
neutrino mass, is given as f*VL which is proportional to Md — Me. The charged-lepton 

mass matrix is written as Me = Cd{Md + k,Mu). We define the complex number ^ as 



^= |Q|(Mrf + /€Mj33/m,. 



(17) 



By definition, |^| ^ 1. The phase of ^ is determined by the charged lepton mass fitting. The 
(3,3) element of Me is then ^e'^'^rrir, where a is a phase of Cd, e*'^ = Cd/\cd\. We denote ^ = ^e*°". 
The pure type II neutrino mass matrix is given as 



M'J oc il-Cd)Md-CuMu 



I 



Cd 



nXdC 



i4>a 



Vubrrib 
Vcbrrib 



\ Vubrrib Vcbrrib {rrii, - im^) / {I - Cd) J 

In the expression, CuTtlc and CuTtIu terms in the diagonal element are neglected. Choosing the 
phase (T to make ^ — 1, the (3,3) entry of the matrix is m^, — m^. Since h-r mass convergence 
implies that and m,- become close to each other as we move to the GUT scale, we can 
assume that m^, — nij. ~ (1 — Cd)Vcbfnbi which implies that the neutrino mass matrix is given as 

/ A2 A A \ 



M,, 



A 
V A 



1 1 
1 1 / 



mo. 



(19) 



where A ~ Vus ~ 0.2. We then obtain bi-large neutrino mixing: both atmospheric and solar 
angle are large. The third neutrino mass, m^g, is about mo. To obtain the large solar mixing, 
determinant of (2-3) block of the neutrino matrix is required to be less than XttLq, and then 
the second neutrino mass, m^^, is about Xuiq. So the mass squared ratio is predicted as 
Amg/Am^ ~ A^. The 1-3 neutrino mixing is also predicted as Ues ~ A. 

Next we consider the case of mixed type II neutrino mass matrix. This case is more compli- 
cated than the pure type II case. From Eqs. (jl3flT^ . the mixed type II seesaw matrix is given 
as 

_y^nmixed ^ _ ^^^^^ ^) ^ 2^^M„ + r^^)' M,(M, - Me)-iM,. (20) 

The first term is just the same as in the pure type II case. When A = 0, we get the type I case. 
The second and third terms might give a hierarchical structure to the (2-3) block of neutrino 
mass matrix, and would spoil large atmospheric mixing. The (3,3) element of the Eq.(j2ni) can 
be written approximately as 

((l_|),l,A),2l^,c(|^)y"., (21) 

where ( = [{Ma — Me)~^]33m^. This element is of the order of m,- or larger in general. However, 
if the (3,3) element can be canceled to the order of (1 + A)(l — CdjVcbrrir, then the bi-large 
mixing neutrino matrix can be realized. Assuming that ^ ~ 1 as in pure type II case, the 
cancellation condition of the second and third terms is q = — (6 + C)/(2 — C). In our numerical 
analysis, the parameter ( is almost real and C — 1- Thus, the cancellation condition is q ~ —7. 
At that time, the phase of Cd is almost vr, and then the phase of ^ must be almost vr from the 
condition ^ ~ 1. The condition can be satisfied if the strange quark has a smaller value of 
mass, since \cd\ ~ m^/ms naively. According to our numerical analysis, ^ is not necessary to be 
1 and Cd can be different values from —7. However, \cd\ has a lower bound, \cd\ ^ 5 to cancel 
the (3,3) element of the neutrino mass matrix in our analysis. Consequently, smaller value of 
strange mass is favored in the mixed type II case. 

In both pure and mixed type II (and of course in the type I case), the neutrino mass 
matrix structure in Eq. p9|) which gives bi-large mixing can be obtained. Then we have a good 
prediction of the 1-3 mixing angle, \Ue-i\ — A. Actually, in the numerical studies in Ref.p21 (type 
II) and Ref.|H] (type I), the prediction of |f/e3| is about 0.17-0.18. Another important prediction 
of the minimal SO(IO) model is the solar mixing angle and the ratio of mass squared differences, 
Amg/Am^. The predictions are almost determined by the (2-3) block of the matrix. The (2-3) 
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block is approximately written as 

where e is a cancellation factor, em-r ^ (M^ — Me)33/(1 — q) in the pure type II case, and em,- is 
approximately Eq.(j2^ divided by 1 — q in the mixed type II case. The |e| should be less than 
Vcb ~ to give rise to a large atmospheric mixing. The condition to obtain a large maximal 
mixing angle is that the determinant of the (2-3) block is canceled, 

\ee"f"msm^ - V^,ml\ ^ OiX)V^,ml (23) 

This condition also provides a small mass differences squared ratio, ArriQ/ Am\ ~ A^. To 
satisfy the condition, one needs ee^'^" ~ +|e|, and |e| can not be very small. Then, maximal 
atmospheric mixing requires a relation — |e|m^ <^ AVcbTrii,. Assuming ^ |e|?7i^, we obtain 
the condition for the bi-maximal neutrino mixing and the small mass squared ratio as 

rris ~ Vcbrrib. (24) 

From the experimental data, the strange/bottom mass ratio at the GUT scale is always smaller 
than Vcb, thus a larger strange mass gives rise to larger mixing angles in the minimal SO(IO) 
model. Furthermore, in the minimal model, the strange quark mass is constrained due to 
electron mass fitting. As we can see in Fig.l, the maximal values of the strange mass depend 
on the KM phase, and the second quadrant KM phase gives larger strange quark mass. As 
a result, the bi-maximal neutrino mixings and the small mass squared ratio favor the second 
quadrant KM phase. This result is in agreement with the numerical studies in Ref. jT^j. In the 
case of mixed type II or type I, smaller values of strange quark mass are favored to obtain the 
cancellation of the (3,3) component of the neutrino mass matrix as we have noted. So in type I 
and mixed type II case, there is a tension between getting large neutrino mixings in general and 
getting their precise values as well as the desired mass square ratio. As a result, the minimal 
SO (10) model with type I or mixed type II predicts that the solar mixing angle is not very 
large and the mass squared ratio is not very small as long as the atmospheric neutrino mixing 
is maximal. We give our numerical data plotting for the solar mixing angle and mass squared 
ratio in the Fig|21 In the mixed type II or type I case, the mixing angle and mass squared ratio 
are bounded for any KM phase, and recent data fitting excludes such bounds more than 2a 
level 0. 
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sin 2e.|2 



Figure 2: The solar mixing angle and ratio of mass squared differences are dot-plotted for 
sin^2^23 > 0.89. The allowed region from recent data fitting lies in the right-bottom corner in 
the figure. 



Next we will see the tan/? bound in the minimal SO(IO) model. From Eq. p5|) . we obtain 



On 



Cd 



cot /?, 



A. 



Cd 



cot p. 



Since and (3u4 are Higgs mixings, we have unitarity constraint la^^dp + \[3u,t 
we have (|an/ad| - l){\(3u/ (3d\ - 1) < 0, and 



mm 



Cd + 3 



1 - Cd 



< tan j3 < max 



Cd 3 



1 - Cd 



(25) 
1. Thus 

(26) 



In the case where ^ ^ 1 which is favored to obtain bi-large structure, (1 — Cd)/c„ ~ mt/mb, 
and then we have lower bounds of tan/5, such as tan/3 ^ mt/mi)\{cd + 3)/(l — Cd)|. In the case 
where Cd ~ —7 which is favored in the mixed type II (and type I) case, the lower bound of tan f3 
is roughly 30. To obtain small tan/5 such as tan/5 ~ 10, Cd ~ —3.5 is needed. Such small \cd\ 
can happen only in the case where the KM phase is close to 180 degree in the minimal SO(IO) 
model. 

In the type I seesaw case, the scale of the right handed Majorana mass is predictable. Using 
the Eqs. fll3llT^ . the light neutrino mass scale, mo in Eq. ()19|) . can be written as 



mo 



'1 + 4/5dt'd 
^1 - Cd) ( ) Vcbrrib- 



(27) 



So using the mass squared magnitude from the atmospheric neutrino data, we can fix the mass 
parameter, mo, and obtain the right handed Majorana mass matrix, Eq. lfT^ . The lightest 
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Majorana mass, which is an important scale for leptogenesis, is about 10^^ GeV. Also we 
obtain the upper bound of vr since \[3d\ is less than 1. 

Now let us summarize the predictions of the minimal SO (10) model. In both pure type 
II and mixed type II, we can obtain the bi-large mixing structure, Eq. (fTI?|) . for the solar and 
the atmospheric oscillation. In the bi-large neutrino matrix, the 1-3 mixing is well predicted 
as If/esl — 0.17 — 0.18. In the mixed type II (and type I) case, larger |q| is favored to obtain 
bi-large structure in the neutrino mass matrix. In order to obtain both bi-maximal neutrino 
mixing and small mass squared ratio, a larger strange quark mass is favored, and this requires 
a smaller |cd|. So, in the mixed type II (and type I) case, the bi-maximal neutrino and the 
small mass squared ratio are not favored, and it is hard to achieve a best fit value of the recent 
combined data analysis. In the pure type II case, the condition of bi-large mixing is just ^ — 1, 
and \cd\ is not constrained. However, due to the electron mass fitting, the strange mass is 
related to the KM phase, and if the KM phase is 60-80 degree, the strange mass can not be 
large and thus it is hard to achieve the bi-maximal mixing and the small mass squared ratio 
in the pure type II case as well. If the KM phase is in the second quadrant and especially for 
the larger value of KM phase, the strange quark mass can be large and the best fit values for 
neutrino oscillation can be achieved. The tan (3 in the mixed type II and type I model is greater 
than 30 roughly and the B — L breaking scale has an upper-bound in the type I model. 

3 Inclusion of Nonrenormalizable terms 

As we have seen in the previous section, the fitting of three charged-lepton masses (especially, 
electron mass) gives relatively small strange quark mass in the case where the KM phase is 
60-80 degree. Then the (2,2) element of the neutrino mass matrix is relatively smaller than 
(2,3) element, and as a result, we get bounds for the solar mixing angle and the mass squared 
ratio. In this section, in order to reproduce the CKM model of CP violation at low energies, 
we employ corrections to the SO (10) model which originate from non-renormalizable terms in 
the superpotential. 

The most simple non-renormalizable term is written including an extra singlet S 

W^y = \-^^{K^i^ii^^H,, + f:^^,^,H,2e) • (28) 

In this case, there is no change from the results of the renormalizable case since the effect of 
these terms is simply to redefine the couplings h and / of Eq.(3). We therefore seek other 
renormalizable terms. 
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We will employ 210 Higgs multiplet to remedy the situation. Note that 210 Higgs (to be 
denoted S) can be used to break the SO(IO) symmetry down to SU{3)c x SU{2)l x SU{2)ji x 
U{1)b-l if we give the (1, 1, 1) and (1, 1, 15) (under G224) submultiplets of 210 vev. One can 
then include the term 

\jfK^^^^^HT. (29) 

Since E has four SO(IO) indices, the combination ifE that can couple to spinors is either a 120 
or 126 representation. We restrict the high scale theory (ultraviolet completion) from which 
the nonrenormalizable term originates in such a way that no 120 term in the effective ifS term 



couples to matter. For instance, we could add a second 126 + 126 pair of field with a high 
mass M ^ Mu which does not develop a VEV with a superpotential given as follows^: 

W' = h'i)i)A' + MA'^' + (30) 

The effective theory below the scale M has a nonrenormalizable term of the type in Eq. (30) 
where only 126 part of Sif field product contributes. If only (1, 1, 1) and (1, 1, 15) submul- 
tiplets of 210 acquire vev, the effective operator that contributes to fermion masses then has 
(2, 2, 15) quantum numbers under (7224- This leads to mass formulae for the quarks and leptons 
of the following type. 

= {h*au + rPu)vu + h'*auVu{^)/Mp, (31) 
Md = ih*ad + r/3d)vd + h'*adVd{^)/Mp, (32) 
Me = {h*ad-3rPd)vd-3h'*adVd{^)/Mp, (33) 
M,^ = {h*au-3rf]u)vu-3h'*auVu{J:)/Mp, (34) 
Mr = rvR, (35) 

where Vu,d are VEVs of MSSM Higgs doublets and vp is VEV of Ap. Note also that h' is not 
a symmetric matrix in general, but we assume that h' is a symmetric matrix like h and / for 
simplicity. 

We can now rewrite the above equations as follows 



Me = M° + AMe, (36) 
^ ^l±^(M,-M°) + M„ + i^AM„ (37) 



'For another nonminimal version with muhiplc 126 at the GUT scale, see Ref.fT^. 
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where = CdMct + c„M^ and AMg = —4:h'*adVd(^)/Mp. In this sumrules, we can break the 
relation of the (2,2) element between charged-lepton and neutrino mass matrices, and thus we 
can increase (2,2) element of seesaw neutrino mass matrix. In the mixed type II (and type I) 
case, we have to increase the (2,2) element while keeping \cd\ to be larger than 5. In fact, we 
find a large solar mixing solution in type I when we switch on the (2,2) element of the AMg. 

In the pure type II case, we do not need the (2,2) element of AMe since a small \C(i\ can give 
rise to maximal mixings. Instead, we need the (1,1) element to cancel and to give the desired 
electron mass and we then find a large solar mixing solution in the case where KM phase is in 
the first quadrant. 

4 Numerical Results 

As already noted, the minimal SO(IO) model is predictive because fermions have only two 



Yukawa couplings, one with 10 and one 126 Higgses. This leads to a sumrule for the leptonic 
mass matrices which determine the the Dirac neutrino and right handed Majorana mass ma- 
trices from the observed neutrino data. In a generic model for neutrino masses, the neutrino 
oscillation data gives us information about the neutrino masses and mixing angles of light neu- 
trino mass matrix, but no information about the Dirac neutrino and the right handed Majorana 
mass matrices separately. On the other hand, a separate information of the Dirac neutrino and 
the right handed Majorana mass matrices are important to extract the predictions of the mod- 
els for lepton flavor violation and leptonic CP violation. In the minimal SO(IO) model however, 
these matrices are completely determined. 

In type I case that we are considering in this paper, we get: 



/ -0.00159 - 0.00014i 0.00067 - 0.0036i 0.017 + 0.015i \ 

0.00369 - 0.00026i 0.0182 + 0.0107i -0.046 - 0.0228i 
V -0.022 + 0.0085i -0.02 - 0.0477i 0.58 + 0.126i y 



(39) 



and 



M, 



R 



( 0.00059 + 0.000048Z -0.00022 + 0.0012z -0.0058 - 0.0051z \ 

-0.00022 + 0.0012i -0.0014 - 0.0063i 0.011 -0.013i 
V -0.0058 - 0.0051i 0.011 - 0.013i -0.037 - 0.0086i y 



X 10^^-°^ GeV. (40) 



The above matrices are calculated for tan /5 = 40 and are in the basis where the charged lepton 
masses are diagonal. The scale 10^^-''^ GeV is the VEV of 126 Higgs which couples to fermions 
and this magnitude of the scale is the maximum possible value given the inputs of the quark 
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masses and the CKM mixings. The h and / in the basis where My^ diagonal are given as follows: 



/ 0.0012 - 0.0002Z 0.00067 - 0.00014Z 0.0011 + 0.0035z \ 
h* = 0.00067 - 0.00014Z -0.0029 - 0.00122z -0.037 + 0.0043z , (41) 
V 0.0011 + 0.0035^ -0.037 + 0.0043i 1.15 + 0.129^ / 

/ -0.0012 - 0.00006i -0.0025 + 0.000062i 0.00088 - 0.0021i \ 
/*= -0.0025 + 0.000062i -0.0097 + 0.00061i 0.021 - O.OOIH . (42) 
V 0.00088 - 0.0021Z 0.021 - O.OOllz -0.0012 - 0.0345z / 

The neutrino mixing angles and mass squared ratio are given as: 

sin^^is = 0.87, sin^ ^23 = 0.92, jt/esl = 0.22, Aml/Am\ = 0.051. (43) 



At tan P = 50, the same matrices are given by: 

/ -0.0017- 0.00042 0.0012 -0.0039i 0.016 + 0.019i \ 

0.004 - 0.00056Z 0.0194 + O.Ollz -0.053 - 0.023z 

-0.023 + 0.01Z -0.025 - 0.052Z 0.660 + 0.13z j 



V 



and 



/ 0.00117 + 0.00026i -0.00076 + 0.0024i -0.010 - 0.012i \ 



R 



-0.00076 + 0.0024Z -0.00099 - 0.0118z 0.0275 - 0.0234z 
-0.010 - 0.012Z 0.0275 - 0.0234Z -0.066 - 0.025i / 



(44) 



X 10 



13.82 



GeV. (45) 



The neutrino mixing angles and mass squared ratio are given as: 



sin^^iz = 0.81, sin2e23 = 0.90, \Ue3\ = 0.22, Am^/Am,^ = 0.06 



(46) 



In the pure type II case, we get 



/ 0.0008 -O.OOlz 



-0.001 - 0.0057Z 0.004 - 0.026Z \ 



0.003 + 0.0013i -0.0056 + 0.00015i -0.053 - 0.008i 
V -0.000039 + O.OOMi -0.012 + 0.00096i 0.627 + 0.1019i / 



and 



R 



/ -0.0016 - 0.000047Z -0.0036 - 0.000029z 0.00077 + 0.0028z \ 

-0.0036 - 0.000029i -0.014 - 0.00044i 0.026 + 0.0008i 
V 0.00077 + 0.00282 0.026 + 0.0008i -0.0158 - 0.0078i / 



(47) 



X 10 



16.24 



GeV. 



(4J 



The above matrices are calculated for tan /5 = 50 and are in the basis where the charged lepton 
masses are diagonal. The scale 10^^'^^ GeV is the scale where SO(IO) gets broken to G2231 which 
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subsequently breaks down to the SM. The neutrino mixing angles and mass squared ratio are 
given as: 



sin^^ia = 0.85, sin^ ^23 = 0.91, \Uez\ = 0.22, Am|/Am^ = 0.084. 
At tan (3 = 40, the same matrices are: 



(49) 



/ -0.00027- 0.00066i -0.0046 + 0.00089i -0.023 - 0.0023z \ 

-0.00073 + 0.00282 0.00074 - 0.0037z 0.007 - 0.052z 
^ -0.0008 - O.OOOlz 0.0013 - 0.0079Z -0.103 + 0.57z / 



(50) 



and 



R 



( -0.00078 + 1.5 ■ 10-S -0.0018 + 0.00040i 0.00039 + 0.0012i \ 

-0.0018 + 0.000040Z -0.007 - 8.7 ■ lO'^z 0.012 + 7.8 ■ 10"^ 
V 0.00039 + 0.0012Z 0.012 + 7.8 ■ lO'^z -0.0081 - 0.00013z / 



X lO^^-^^ GeV. 



(51) 



The neutrino mixing angles and mass squared ratio are given as: 



sin^^ia = 0.85, sin^ ^23 = 0.92, \Uez\ = 0.20, Am|/Am^ = 0.058 



(52) 



Using these fits, we will calculate the lepton flavor violating processes and the amount of 
leptonic CP violation in type I and pure type II scenarios in the mSUGRA model in sec. 6. 

Even though all the examples we have given have large Ues, its value can be much smaller 
due to the presence of new parameters in the higher dimensional term. 



5 Gauge Coupling Unification 

From the results of the previous section (Eq. (40) and (45)), we see that for the type I model to 
be successful in predicting neutrino masses, the intermediate scale (vr) must be below the GUT 
scale and around 10^^ GeV for our examples. We assume the resulting symmetry breaking chain 
is of type, ^O(IO) ^ SU{3)c x SU{2)l x SU{2)r x U{1)b^l ^ SU{3), x SU{2)l x U{1)y. 
This is very different from the gauge coupling uniflcation scenario in MSSM. It is therefore 
important check if the type I models are compatible with gauge coupling unification. 

In our model, in the scale region Msusy < < vr, the spectrum is that of familiar MSSM. 
Above the vr scale, the symmetry group expands to SU (3)c x SU (2) l x SU {2)r x U{1) b~l- In 
addition to the new gauge bosons associated with this new symmetry, the new matter and Higgs 
that contribute are the following: three RH neutrinos as part of the SU {2)fi lepton doublet; one 
bidoublet from which the MSSM doublets emerged. In addition we need either the B — L = 2 
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Figure 3: The couplings l/oj are plotted as a function of scale. 

triplet pair (A^ + © A^^ + Ar) and x[(3, 3, 1, -2/3) + (3, 3, 1, 2/3)] © t/[(3, 1, 3, -2/3) + 
(3,3, 1,2/3) to have successful gauge coupling unification. If we use the Higgs spectrum from 
Ref. [T^ in the context of our model, we find that x ®y (from 210) and A^ + © A/j + A/j 
are at about the vr scale. Using these fields, the Fig. El shows the coupling unification in this 
model and the gauge unification happens at about 10^^'^ GeV. We have chosen the x, y and A 
Higgs masses to be 6 times the vr scale to get this value. 

First of all it is gratifying that the type I model is compatible with gauge coupling unification 
with the desired value of vr scale. Furthermore the lowering of the ultimate unification scale 
implies that the gauge exchange contributions to the proton decay in this model lead to proton 
lifetimes much lower than the MSSM unifying without an intermediate scale and will be around 
the current lower limit 5 x 10^^ yrs when we include the threshold corrections |T7j. The profile 
of the proton decay modes in the type I case will be substantially different from the pure type 
II case recently discussed fH] and is presently under investigation. 

6 Lepton Flavor Violating Processes and Lepton Edm 
6.1 Lepton Flavor Violation 

We now discuss the lepton flavor violating processes e.g. /i 67, r — > /i7 etc. The operator 
for li Ij + 1 is: 



Ci^^i^^ = - — Ij cr^^qu {aiPi + a^P/?) k- + h.c. (53) 
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where Pl,r = {IT 75)/2 and ct'^^ 



I [7^, 7^^]- The decay width for /j ^ + 7 can be written 



as: 



647r 



) 



(54) 



Then the branching ratio is obtained by multiplying this decay width with the life time of the U 
lepton. The supersymmetric contributions include the neutralino and chargino diagrams [TH] . 

We work in the basis where the charged lepton masses are diagonal at the highest scale 
of the theory. We first start with pure type II seesaw model. We assume that vr (the scale 
where G2231 breaks down to SM) is at the GUT scale. So we have just MSSM and right handed 
neutrinos below the GUT scale. The right handed masses have hierarchies and therefore get 
decoupled at different scales. The flavor-violating pieces present in Yy induces flavor violations 
into the charged lepton couplings and into the soft SUSY breaking masses e.g. terms etc. 
through the following RGEs: 



We use the mSUGRA universal boundary conditions at the GUT scale to draw the figures. 
In Figs, mini we show the BR[/i 67] and BR[r /i7] as a function of mi/2 for different 
values of Aq. The lightest neutralino is the dark matter candidate in this model and we satisfy 
the 2(T range of the recent relic density constraint ^cdm = 0.1126]'^o;oQ9j2n] in the parameter 
space. When we satisfy the relic density constraint, the mo gets determined. For example, mg 
varies between 60-100 GeV for = GeV line in the graph. The Figs also show that the 
larger tan (3 has larger lepton flavor violation. 

In type I model, the upper bound of Vr, VEV of 126 which couples to fermions, is determined 
and this value is not the GUT scale. For example, vr can be < 10^^ GeV. Then, (72231 symmetry 
can be maintained between the GUT and the vr scale if other Higgs fields do not break B — L 
symmetry. This feature induces larger lepton violation since the right handed neutrino is a 
part of the doublet which has right handed electrons. The right slepton masses get new flavor 
violating contributions through the flavor violating pieces present in Y^. The new contributions 
to the RGEs: 



dYjdt 




(55) 





dYjdt 




(56) 
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Figure 4: The BR[/i 67] is plotted as a function of mi/2 for different values Aq and tan /3 = 10, 
40 and 50 in pure type II. 



dm\j^/dt 



167r2 



The new effects arising from the RGEs make the lepton flavor violations to be larger in this 
case and is depicted in Fig. El In Fig. IHl we plot t fi'j and find that the Br can be as large 
as 10~* which can be explored in the near future. 

6.2 Electric Dipole Moment of Electron 

The EDM, dj for fermion fi appears in the effective Lagrangian as: 



(57) 



We have contributions from the chargino and neutralino diagrams to (i/[2I]- The electron EDM 
is plotted in Fig. [7|for the type II and in Fig. |Hlfor the type I. We again use the same SUSY 
parameter space as in the case of lepton flavor violation. We find that the maximum value of 
EDM is ~ 10~^^ ecm for the type II. For the type I, the EDM is large and is around 10~^^ ecm 
and a large region of parameter space can be explored very soon. The muon EDM is shown 
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67] is plotted as a function of mi/2 for different values Aq and tan /3 = 40 
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Figure 6: The BR[r — > /i7] is plotted as a function of mi/2 for tan/? = 40 in pure type II and 
type I. 
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Figure 7: The electron EDM is plotted as a function of mi/2 for different values Aq and 
tan /5 = 40 and 50 in pure type II. 



in Fig. ini and the maximum value shown is about 10~^^ ecm and the scaling is broken in this 
model. 

The sin 2/5 calculated from the —>■ (pKg mode in this model is 0.67, which is not different 
from the SM prediction. 



7 Conclusion 

In summary, we have revisited the minimal renormalizable SO(IO) models with a 126 Higgs 
multiplet that has recently been shown to predict neutrino mixings in agreement with ob- 
servations, with the primary goal of reconciling CKM CP violation with successful neutrino 
predictions. We consider the most general type II and type I seesaw formula for neutrino 
masses that includes the effect of the right handed neutrino mass matrix. We find that in these 
models the basic ingredients that went into understanding maximal neutrino mixings i.e. the 
relation M^, = c{Md — Me), can be recovered in certain limits. However they do not help to 
keep the CKM phase in the first quadrant and in the type I case are incompatible with detailed 
neutrino data. We remedy this by including a specific class of nonrenormalizable terms in the 
Yukawa superpotential, that follows from a simple high scale theory. This specific set of non- 
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Figure 8: The electron EDM is plotted as a function of m.1/2 for different values Aq 
tan P = 40 and 50 in type I. 
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Figure 9: The muon EDM is plotted as a function of mi/2 is plotted as a function of mi/2 
40 in pure type II and type 1. 
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renormalizable terms also reconcile pure type II models with the CKM model of CP violation. 
We find that the type I model requires an intermediate vr scale (i.e. vr <^ Mjj) which is lower 
than the standard GUT scale, however it is compatible with gauge coupling unification which 
happens around 10^^'^ GeV. The proton life time is around the current experimental limit. 
We also find that the mixed type II and type I models require a value of the supersymmetry 
parameter tan/3 larger than 30 to be compatible with present neutrino data. We then study 
the phenomenological implications of the type I, the mixed type II and the pure type II models 
for neutrino mixings, lepton fiavor violation and lepton edms. We find that these predictions 
provide a new way to test these models. 
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